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We analyze the effects of a random magnetic potential in a microfabricated waveguide for ultra- 
cold atoms. We find that the shape and position fluctuations of a current carrying wire induce a 
strong Gaussian correlated random potential with a lengthscale set by the atom-wire separation. The 
theory is used to explain quantitatively the observed fragmentation of the Bose-Einstein condensates 
in atomic waveguides. Furthermore, we show that nonlinear dynamics can be used to provide 
important insights into the nature of a quantum phase transition from the superfluid to the insulating 
Bose glass phase may be reached and detected under the realistic experimental conditions. 
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Several groups have recently reported realizations of 
quasi one-dimensional (QID) Bose-Einstein condensates 
(BEC) in the magnetic waveguides for atoms created by 
the microfabricated wires (magnetic traps on microchips) 
[1234]. A surprising feature is the presence of fragmen- 
tation: a random modulation of the atomic density in 
the axial direction . For smaller atom- wire separa- 

tion the fragmentation is enhanced and its characteristic 
lengthscale decreases. Determination of the origin and 
properties of the observed density modulation is essen- 
tial for understanding atomic condensates on microchips 
and for developing atom-optical microfabricated devices. 
Several observations suggest that the wire shape and po- 
sition fluctuations are the primary origin of the fragmen- 
tation l^l^lll, however, a detailed theoretical understand- 
ing of these phenomena is still lacking . 

In this Letter we use a first-principle microscopic cal- 
culation to demonstrate that a small meandering of the 
wire leads to a strong random potential in a magnetic 
microtrap Q . We show that even when the wire fluctua- 
tions have no intrinsic lengthscale, the disorder magnetic 
potential in the waveguide is a Gaussian correlated ran- 
dom potential with a correlation function that is strongly 
peaked at a finite wavevector and vanishes for small and 
large wavevectors. The characteristic length scale of this 
disorder correlation function is set by the atom- wire dis- 
tance, d. This theoretical model is used to quantita- 
tively explain the experimental results. We then study 
the properties of the QID atomic BEC in the regime of 
strong fragmentation, which may be realized in micro- 
traps by bringing atoms closer to the chip surface and/or 
changing the atomic density. Some of the most strik- 
ing manifestations of this strong disorder limit appear in 
the non-linear dynamics of the fragmented condensate, 
including a chaos and a self-trapping of the excitations. 
We argue that the quantum phase transition between 
the superfluid (SF) and the insulating Bose glass (BG) 
phases can be achieved and detected under realistic ex- 
perimental conditions. 

Neutral atoms with a magnetic dipolc moment, jia, 
anti-aligned with respect to the magnetic field, B, ex- 
perience a potential ^(r) = /iJB(r)|, so they can be 
confined near a field minimum [J. A typical microtrap 
setup is shown in Fig. ^ The waveguide is located at 
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FIG. 1: Typical setup of a magnetic microtrap: A current Jo 
flowing in a microfabricated copper conductor and a perpen- 
dicular bias field, B±, form a magnetic waveguide for atoms. 
An offset field -B|| applied parallel to the wire reduces loss 
processes near the magnetic field minimum. 



a distance d = -^s. from the wire center, where the az- 

cB± ' 

imuthal magnetic field created by the wire cancels out 
the transverse bias field, B±. The longitudinal confine- 
ment can be provided by the additional magnetic field 
gradients applied parallel to the wire. 

We now discuss a random magnetic potential along 
the center of the waveguide caused by the shape and 
position fluctuations of the current carrying wire. Let 
fL/B.{z) be the deviations of the wire's left/right bound- 
aries from their ideal positions at y = ^Wo/2. Changes 
in the wire height provide a much smaller contribution 
to the fluctuating magnetic field, so we will neglect them 
throughout this paper. In a steady state the devia- 
tion of the current density, 63{y,z), from its average 
value, Jo, satisfies the charge conservation condition, 
dSJy/dy + dSJz/dz = 0. If we simultaneously assume 
a constant conductivity throughout the wire, the current 
is also vorticity free, d5Jz/dy — dSJy/dz — 0. Boundary 
conditions for 53{y, z) follow from the requirement that 
no current flows out of the wire. Assuming small fluc- 
tuations of the wire (|/l/_r(2:)| < Wq and \^^^\ < 1), 
we have 5Jy{'^Wo/2, z) / Jo — dfL/i^{z)/dz. It is conve- 
nient to introduce an auxiliary scalar potential Gj{y, z), 
such that SJy = JodGj/dz and 6Jz — —JodGj/dy. The 
function Gj{y, z) satisfies the Laplace equation in the in- 
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terior of the wire and has Dirichlet boundary conditions 
at y = =fWo/2- Using separation of variables, we find 



Gjiy,z) = 



'■dk 



IT J sinh(fcVFo) 

sinh(A:2/)cosh(fciyo/2)4 



cosh(fc?/) sinh(fcW^o/2)/fe 



(1) 



where = 5 dze''-^''{fR{z) ± /^(z)) are the Fourier 
components of the wire position and width fluctuations. 
The fluctuating part of the magnetic field can be found 
from the Biot-Savarat law: 
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To obtain Eq. lI2Jl we performed the integration by parts, 
which led to a cancellation of the contributions from the 
y G {-^ + fdz),-^) and y G + Mz)) 

regions. The leading order random potential comes 
from the z-component of the magnetic field, 5U{t) — 
±fj,aSBz{d, 0, z), if the condensate width is much smaller 
than d (the ± sign is the orientation of the offset field, 
_B||). From Eqs. Q-© one can see that only the wire 
position fluctuations, /+, contribute to 6U{r) at this or- 
der. 

To elucidate the nature of the resulting random po- 
tential we compute the correlation function Aj, — 
J (ize'''(^~^')((5[/(z)JC/(z'))dis, where (■ ■ •)dis is disorder 
ensemble average. Assuming Hq <C d (this condition is 
typically satisfied in experiments), we obtain 



A,- = 



2Io^la\ {kdY 
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where F,. = / dz e-*'=(^-^') (/+(z)/+(z'))dis, and 

2sinh(x/2) ^ (-1)" 
a;smh(a;)Ai(y) n!(2y)" 

X [7(2n + 1, a;/2) - 7(2n + 1, -x/2)] . (3) 

Kn{x) is the Bessel function of the second kind, and 
j{n,x) = dx' x''^~^e~^ is the incomplete Gamma 
function. The strength of the disorder potential can 



be defined as u1 
d ^ Wo we obtain 
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FIG. 2; (a) Correlation function of the random magnetic 
potential in a microtrap, Ak, assuming white noise fluctu- 
ations of the wire position (i.e. Fk — constant). Insert: 
Afc when the wire fluctuations have an intrinsic length scale: 
Fk (X (e-('=-*i)''''? +e-(fe+'=i)''''?). We take 27r/A;i = 200^m, 
rii = lOO^m, and (/+(z)/+(z)>dis = (0.1 fimf. Solid to 
dash-dotted lines correspond to d = 50, 100, 150, and 200 /im 
respectively, (b) Condensate density profiles with parameters 
chosen to be close to the values used in the experiments of Ref. 
[1^. Dotted lines are the results without random potential. 



Fixing current Iq we find Ug ^ d^^^^, which is close to 
the experimentally observed, Ug ~ d^^-^ ^.We can also 
estimate the magnitude of the random potential: If 6f ^ 
0.1 /im is the scale of the wire position fluctuation and 
^ ~ 100 /zm is the appropriate correlation length, then 
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kHz for ^ 1 G and d - 100 fim. 



This value is comparable to a typical value of the EEC 
chemical potential and explains the strong fragmentation 
effects observed in magnetic microtraps |33ll4lla| . 

In Fig. |2Ia) we show Afc computed for the white 
noise fluctuations of the wire center position It peaks 
at ko « 1.3/d and vanish in the fc = and k — 00, 
showing a characteristic lengthscale, d. Long wavelength 
fluctuations are suppressed because a uniform shift of the 
random potential does not lead to a parallel (z) compo- 
nent of the magnetic field. In the insert of Fig. HJa), we 
demonstrate how the random magnetic potential is mod- 
ified when the wire fluctuations have an intrinsic length- 
scale 1/ki (see caption) . When the wire fluctuations are 
short ranged (kid 3> 1), the potential fluctuations re- 
main peaked close to ko ^ 1.3/d, and the lengthscale of 
the random potential is set by the atom- wire separation. 
When the wire fluctuations are longer ranged (kid < 1), 
Afc is then peaked at k ki, and the random potential 
tracks the wire fluctuations. A linear relation between 
the condensate height and the lengthscale of the frag- 
mentation was observed in Ref. 3 4]. In Fig. Olb) we 
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FIG. 3: (a); Schematic dynamical phase diagram of frag- 
mented quasi one-dimensional BEG in relation to the shaking 
experiments. The dashed line indicates the minimum dis- 
placement to generate an observable density fluctuations (i.e. 
C > Cmin). Three types of the system response can be iden- 
tified from the power spectra of the oscillations: (b) dipole 
mode regime, (c) chaos, and (d) self-trapping of the excita- 
tions. Power spectra are obtained from simulations of a three 
well system with {Nj ) = lO'^ . 



show the calculated condensate density profiles for 10^ 
sodium atoms in the Thomas-Fermi approximation using 
the same wire fluctuations taken from the distribution of 
the insert of Fig. [2Ia). The fragmentation appears for 
d < 100 /im (c.f. Ref. 3]). We note that the profiles 
obtained by a white noise potential (i.e. =const.) 
do not have a characteristic length scale, and are very 
different from the above results and experiments. 

We next consider the new phenomena that emerge in 
the strong fragmentation regime, where the random po- 
tential is appreciably larger than the chemical potential 
of the atoms so that the tunneling rate for the atoms 
between the neighboring wells is much smaller than the 
local confinement frequency. The zero temperature con- 
densate wavefunction then can be approximated by the 
tight-binding model $(z) = J2j y^^'ji^) ^ 

where (j^ji^) is the localized single particle wavefunc- 
tion (here obtained by solving the onsite Gross-Pitaevskii 
equation), Nj is the number of particles, and Sj is the 
phase of the mini-condensate j. Equations of motion of 
the mini-condensates (i.e. fragments) are given by 



j+a 




0, 



(5) 
(6) 



where Uj and Nj are the "charging energy" and the 
equilibrium number of atoms in the well j; Kj ji = 
Kj_jt ^yNjNj' is the "Josephson energy" and Kjj> is a 
single particle tunneling rate between the wells j and j' . 
Eq. © corresponds to the Josephson relation between 
the local chemical potential S^j = Uj{Nj — N^) and the 
phase Sj. Eq. ^ describes the charge conservation. 
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FIG. 4: The dotted line shows the probability P{Q < 1), 
which is close to zero if the system is in the superfluid regime 
and is close to one if in the insulating Bose glass phase. The 
solid line shows the probability that one can observe fluctua- 
tions in the number of atoms between the neighboring wells 
with LOj > tJmin and > ^min (sce text). When this probabil- 
ity approaches zero, the system will appear self-trapped in the 
shaking experiments. A'^avo is the average number of atoms in 
a single well (mini-condensate). Disorder strength for (a) is 
the same as used in Fig. Hlb) for d = 100 /im (denoted to be 
Mj), while it is five times stronger in (b). 



The key properties of such fragmented BEG can be re- 
vealed by exciting BEG in the "shaking" experiments. In 
these experiments one quickly displaces the global con- 
fining potential along the waveguide and observes the 
ensuing motion of the condensate. From the numerical 
analysis of Eqs. lO and ©, we identified three types 
of response of the system to the initial displacement, D, 
(see Fig. Ela)) 9]: (i) For small D only the dipole mode 
is excited, (ii) As D increases, the dipole mode becomes 
unstable and many other modes are generated, leading 
to a chaotic behavior, (iii) When D gets even larger, the 
system is "trapped" in small oscillations around some 
configuration with strong imbalance of the atom density. 
The number of atoms in individual wells is not oscillat- 
ing around the equilibrium values N^, but around some 
non-uniform density distribution. We call this regime 
self-trapping of the excitations ■ 

The power spectra of the oscillations provide a direct 
way to identify different types of the dynamical behavior. 
In the regime (i) condensate oscillates at the frequency of 
the dipole mode and its higher harmonics (see Fig|3fb)). 
In the regime (ii) we have a broad distribution of frequen- 
cies (see Fig|2Ic)). In the limit of self-trapping, for every 
point in the system we have several dominant frequencies 
and their harmonics (see Fig|2|^d)). The transition of the 
system from a dipole mode into a chaotic regime origi- 
nates from the non-linearity of Eqs. (0-® and is similar 
to the modulation instability in optical lattices It is 
commonly believed that for strong disorder the superflu- 
idity are terminated by strong quantum fluctuations and 
that the system should go into the insulating Bose glass 
phase (BG) Therefore, the coherent dipole fluctu- 
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ations should disappear at the SF/BG transition point, 
and the self-trapping response should become dominant 
in the BG phase due to the localized nature of the insu- 
lating state. The qualitative phase diagram is shown in 
Fig. Ha). 

More quantitatively, transition into an insulating state 
can be estimated by the ratio of the "Josephson energy" 
between the neighboring wells to the average charging 

energy, = ^^^^^^^-^^'-^^^ . The SF to the Mott insula- 
tor transition in a periodic potential takes place when 
Q ~ lin the meanfield theory. For a correlated random 
potential we expect that the SF to the BG phase transi- 
tion takes place when the probability to have Q < 1 in 
a junction, P(Q < 1), is of the order of one. In Fig. 0] 
we show V{Q < 1) as a function of the atom density for 
two different strengths of the disorder potential (dotted 
lines). Taking the parameters corresponding to the ex- 
periments of Ref. |3 with d = 100 /im, we estimate that 
a system of atoms with A^avc < 10^ per fragment should 
be in the BG phase. 

Naturally, the shaking experiments can be used to 
identify the SF/BG transition. In practice, a finite life- 
time of a condensate limits experimentally measurable 
Josephson frequencies ujj to be ujj > ujmin- When a typ- 
ical uj becomes smaller than Wmin, oscillations may no 
longer be detected. Another limiting factor for the exper- 
iments comes from the finite density contrast, C — 
required for resolving oscillations. When typical ( be- 
come smaller than (^rnin (for a fixed D), the system will 
always appear trapped. The solid curve in Fig. ^ shows 
the probability of observing a dipole mode or chaotic dy- 
namics for a given disorder strength with the minimum 
displacement for uj^i^ = 2tt ■ IHz and Cmin = 0.1. The 
observation of the self-trapping transition as a function 



of the disorder strength at a small D should provide a 
reasonable estimate of the SF-BG transition point 9||, as 
illustrated by the dashed line in Fig. Ola). We expect 
such dynamical properties and a sharp crossover to the 
insulating regime are still valid, at least qualitatively, in 
the temperature regime of current experiments. 

In the shaking experiments discussed so far the random 
potential was kept fixed and the confining potential was 
displaced. It would also be interesting to keep the con- 
fining potential fixed and oscillate the random magnetic 
potential at some small frequency. This would be a direct 
analogue of the Andronikashvilli experiments for super- 
fiuid ^He in a disordered medium |l2l |. which provide a 
direct measure of the superfluid condensate fraction. 

In summary, we performed microscopic calculations of 
a random magnetic potential for the atomic BEG in a 
microtrap. We showed that a small meandering of the 
wire is sufficient to explain the experimentally observed 
fragmentation of the condensates ||6j. The response to 
the trap shaking can be used to study the superfluid 
properties of the strongly fragmented condensates and 
to identify the quantum phase transition from superfluid 
to the insulating Bose glass phase. More generally, our 
work demonstrates that atoms in the magnetic micro- 
trap is a very promising system for controlled study of 
one-dimensional disorder problems. 
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